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I. INTRODUCTION 

Mesoscopic physics concerns with the properties of small systems with sizes in the range from of a few nanometers to 
micrometers, at low temperatures typically below 1 K. The constant progress in nano-fabrication techniques allows for 
a controlled realization of these structures and a consequent increasing interest in this physics. Hybrid superconducting 
' systems are mesoscopic devices constituted by normal and superconducting parts in electrical contact. Characteristic 
5r ! °f superconductivity is the macroscopic phase coherence of the order parameter and the supercurrent flow. On the 
one hand, superconductivity adds new degrees of freedom and makes the physics of mesoscopic systems richer. On 
the other hand, superconducting properties arc deeply influenced by mesoscopic effects as it is the case of charging 
effects in small superconducting junctions. 
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This brief introduction makes evident that the field of hybrid mesoscopic superconducting systems is vast and 
diversified. It is our intention to give a general overview of the field, pointing to the relevant references for a more 
detailed analysis, and to keep the presentation as simple as possible. The paper is divided in five sections (including 
this introduction) which are almost independent among themselves. We begin, in Section II, by introducing the 
basic properties of hybrid systems whose electrical transport is dominated by the conversion of normal current into 
supercurrent which takes place at the interface. We give, for educational purposes, an introduction to the technique of 
the Bogoliubov-de Gennes equations combined to the scattering approach for the study of transport in hybrid systems. 
The complementary approach using quasiclassical Green functions is then briefly introduced, however we refer to the 
existing literature for a more detailed presentation. In Section III we introduce the heat-transport properties of 
hybrid systems in order to discuss the effectiveness of such systems as electronic coolers. In Section IV we provide a 
presentation on the superconducting properties in small systems where charging effects are dominant. 

Basics in superconductivity and Josephson physics can be found in the books by Tinkhaml and by Barone and 
Paterno^. Since many ideas discussed here belong to the field of more traditional mesoscopic physics we refer for 
these topics to the books by Beenakker and van Houter*^ and Imry^ and to the conference proceedings^. Excellent 
introductions to single charge tunneling can be found in Refs.— We also refer to the recent reviews on hybrid 
structures^ and ultra-small superconductors^ for a more complete description of these topics and a more complete 
list of references. 



A fascinating aspect of mesoscopics is the possibility of fabricating hybrid nanostructures formed from combinations 
of normal conductors (N) and superconductors (S). The interest for these systems stems from the fact that the nature of 
charge transport is dramatically different in the normal and in the superconducting constituents of the heterostructure, 
giving rise to a rich variety of effects. While the charge current in N is carried by quasi-particle excitations (electrons 
in a metal), at sufficiently low energy (e.g. temperature, applied voltage, frequency,...) the current in S is carried 
by the superconducting condensate, which is a many body ground state property of the system, and flows without 
dissipation. 

The matching between these two different charge transport mechanisms can be well understood for phase-coherent 
structures in terms of the so-called Andreev reflection process^. This can be introduced starting from the microscopic 
BCS theory of superconductivity^. In the presence of boundaries and non-uniformities it is convenient to perform 
an approximate Fourier transform of the BCS Hamiltonian. Making use of the Hartree approximation, one arrives to 
an effective single-particle Hamiltonian which leads to the Bogoliubov-de Gennes equation^ 



for the coherence factors u and v, where the energy E is measured from the condensate chemical potential. This is 
a matrix eigenvalue equation which contains the single-particle Hamiltonian Hq describing electrons in the supercon- 
ductor in the absence of the attractive potential, and the time- reversal of such a Hamiltonian (— Hq), which describes 
the hole degree of freedom. The off-diagonal term A, which introduces a coupling between electrons and holes, is given 
by the product of the point-like electron-electron attractive potential and the pairing amplitude describing supercon- 
ducting correlations among electrons. The pairing amplitude, proportional to the anomalous average (ipip) can be 
understood as being the wave function of the Cooper pairs. Since A, which is the superconducting order parameter, 
is determined by the coherence factors, it should be calculated self-consistently. In the absence of superconductivity 
(A = 0) the system is described in terms of decoupled electrons and holes, the latter being characterized by having 
group velocity and wave-vector in opposite directions, as well as opposite charge and spin with respect to electrons. 
When an electronic attractive interaction is present (A / 0), electron and hole degrees of freedom get mixed forming 
particle-like and hole-like quasi-particles. As a result, a gap of amplitude |A| opens up in the energy spectrum (see 
Fig. []J,) forbidding quasi-particles excitations with energy inside the gap to propagate in the superconductor. 

Let us now consider a piece of normal metal in contact with a superconductor and assume that electrons of a given 
energy E < |A| are injected from the normal part toward the NS interface. As explained above, such electrons cannot 
propagate through the superconductor, in fact they are allowed to penetrate as an evanescent wave up to a depth 
of the order of the BCS superconducting coherence length £. They can however undergo two different processes. 
In the first one, normal reflection, electrons are reflected back into the normal slab and do not contribute to the 
charge current (sec Fig. [TJd). The second one is made possible by the electron- hole coupling term of the Bogoliubov- 
de Gennes Hamiltonian and consists in the coherent evolution of the incoming electron into a retro-reflected hole 



II. HYBRID SYSTEMS 



A. Andreev reflection and proximity effect 




(1) 
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(see Fig. QJ). This constitutes the Andreev reflection process, which is indeed responsible for charge transport, 
and corresponds to the transfer of a pair of electrons from the normal side of the interface to the superconducting 
condensate. Noteworthy, each Andreev reflection process contributes of 2e to the charge current. It is important to 
notice that Andreev reflection does not contribute to the transport of energy, since no excitations are transferred into 
the superconductor. While providing good charge conductivity, NS interfaces prove to be good thermal insulators. 

Andreev reflection also accounts for the influence that the superconductor exerts on the normal conductor, known 
as proximity effect. In much the same way as the quasi-particlc excitations penetrate the superconductor, the super- 
conductor pairing amplitude also leaks out into the normal side of NS contact over a large distance, of the order of 
the superconducting coherence length, even in the absence of an attractive interaction. As a result, partial supercon- 
ducting properties are induced into the normal conductor. 

In the next Section (QTBj we will give a derivation (following RefJ^) of the Bogoliubov-de Gennes equation, while 
in Section Hi CI we shall solve such an equation for different systems. In the following sections we shall describe the two 
most widely used theoretical approaches applied to study hybrid systems, namely the scattering formalism (Section 
IIIDj) . which we will analyze in details, and the quasi-classical Green's function theory (Section fllFj) . In Section Hi Gl 
we shall review the most important results attained in the field. 



B. Derivation of the Bogoliubov-de Gennes equation 

As mentioned above, the Bogoliubov-de Gennes equation is derived from the BCS theory. In this section we 
show how to derive such an equation in the general case where an exchange field is present. This will be useful for 
studying hybrid structures in which ferromagnets are present too. It is convenient to make use of the Bogoliubov's 
self-consistent field metho d ' 15 . The starting point is the following Hamiltonian in the second quantization form: 



H 



dr U\(r)^\(f))H{7) 



V 



^2 / d ? ^(r)^,(rO^ CT '(r)<Mr) , 



(2) 



where the second term accounts for the electron-electron coupling, which is assumed to be a point-like, two-particle 
interaction (introduced for the first time by Gorkov^) and independent of spin, therefore characterized by a single 
coefficient V. In Eq. ©, 



H(r) 



1 

2m 



c / 



^ he 
= 2mc 



(3) 



is the single-particle Hamiltonian in the presence of an exchange field (second term). A(r) is the vector potential, 
Vo(r) is the normal potential, (i is the chemical potential, M is the internal mean field due to the exchange interaction 
and Uo(r) is the periodic potential due to the ions in the crystal, c is the speed of light, 1 is the 2x2 unity matrix and 
a is the vector of Pauli matrices. Before proceeding, we want to remark that ferromagnetism is introduced according 
to the Stoner model (see Refs . 17 ' 18 ' 19 ' 20 ). In this model a molecular field, due to the permanent magnetization present 
in a fcrromagnct, replaces the external magnetic field in the Zceman energy, and it is added to the single-particle 
energy. Such a molecular field can be obtained calculating the exchange interaction in a system of spins using the 
mean field approximation (see Ref.—). Such a molecular field can be treated as an adjustable parameter to fit band 
structure calculations. 

The operators '(/v(r) and 4>],{f) are, respectively, destruction and creation field operators for a particle of spin a at 
the point r, which obey the Fermi commutation relations 



^(r)ip*>(r') +ipa>(r')iP, y (f) = 
ipU^a'ir') + i><r>(r')i>t(r) = 6 a *>6(r- ?) 



(4) 
(5) 



We now seek an approximate solution of the many body problem in terms of an effective single-particle Hamiltonian 
H e ff. According to the Hartree approximation (see Ref.—), such a Hamiltonian describes a particle which moves in 
a mean-field potential produced by all other particles. Using the variational principle, it can be shown that 



eff = f dr(^\{r)^\(r)) 



H(r) 



W{r) 
U l {f) 



dr 



^{(r)^{(r) A(f) +i' l (f)^(r) A* (f) 



(6) 
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where W(r), U^{r) and A(r) are effective potentials defined by 

W T (0 = -V(^I(f)^ i (f)) , 
U l {r) = -y(^](f)^ T (f)> 



and 



A(f) = -F(^(r)^ T (f)) = y(^ T (r)^(f)) 



(7) 

(8) 

(9) 



W{r) and U^(r) are single particle potentials (standard Hartree result for point-like interactions) and A(r) is referred 
to as the pairing potential. 

The key feature of the mean-field theory of superconductivity is that we admit non-vanishing expectation values 
for the spin-paired operators ipi(f)ipi(f) and ipi(f)ipi(f). If we consider H e ff in ([6]), we see that, on the one hand, 
the terms in and destroy and create one electron and therefore conserves the number of particles. On the other 
hand, the terms in A increases or decreases the number of particles by 2. This is not a problem, since H e ff operates 
on the BCS wave function which is not an cigenfunction of the number operator. Applying the following unitary 
transformation 



(10) 



H e f f can be diagonalizcd in such a way that 



H eff - E g 



(11) 



where E g is the superconducting ground state energy, E n is the energy of the ra-th quasi-particle excitation and % 
{f)Xa) is the quasi-particle destruction (creation) operator satisfying Fermi commutation relations: 



(12) 
(13) 



The Bogoliubov-de Gennes equation is an eigenvalue equation which allows one to calculate the energy spectrum 
of the quasi-particles and the coherence factors u^(r) and v°(r): 



f(H l -h z ) 



V A*(r) 











(H l + h z ) A(f) 
A*(f) (-H 1 +h z ) 
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(14) 



where H (f) = H +U a (f), with 



H n = 



2m 



Vb(r)-/x 



(15) 



he 

h, = —M z , (16) 

having set A = and M x = M y = 0. 

We shall refer to the matrix Hamiltonian in (fl"4|) as the Bogoliubov-de Gennes Hamiltonian Hbg- If we identify 
the top-left 2x2 block of such a Hamiltonian with the particle Hamiltonian H p and the bottom- right 2x2 block with 
the hole Hamiltonian Hh, we see that the equality Hh = —H* holds. 

So far, we have diagonalizcd H e ff using the transformation (fTQ|) . which was proposed independently by Bogoliubov^ 
and Valatin^i, known as Bogoliubov-Valatin transformation. H e f f is now written in terms of quasi-particle operators 
7, 7 t , which greatly simplifies calculations. Although the Bogoliubov-de Gennes equation (fT4|) does not determine the 
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potentials W (r), U^(r) and A(r), they can be fixed self-consistently. By replacing the field operators in the definition 
of the above potentials ([7]), ©, © and using the Bogoliubov-Valatin transformation (JTQJ) we obtain 

^(r) = -y^[|4(f)| 2 / n + k(f)| 2 (l-/„)] , (17) 

n 



U l(f) = Hmffn + l^(r)| 2 (l - /«)] (18) 

and 

n 

where 

/n = • (20) 

e K s T +1 

Equations (jTTj> , fl8]) and (p~9j) constitute the self-consistent equations for W(r), U^(f) and A(r). These ensure that 
if we solve the Bogoliubov-de Gennes equation (fT4|) . then the value of W, and A which is calculated from the 
solutions of (fTTjl . (O and (fT9)) is equal to the initial value of W , and A. 

There is an important distinction between U^{r) and lA^(r), on one side, and A(r) on the other. The Hartree 
potentials (f) and U^[f) come from a sum involving all states below the Fermi level and hence are nearly temperature 
independent and can be approximated by the Hartree potentials calculated in the normal state. However, the pair 
potential A(r) is a sum of terms of the form u°{r) *(r) which are non-zero only in the neighborhood of the Fermi 
surface. For this reason A(r) is a strong function of temperature. 



C. Solutions of the Bogoliubov-de Gennes equation 

In this section we discuss the solutions of the Bogoliubov-de Gennes equation (|14[) in heterostructures containing 
normal metals, ferromagnets and superconductors. In the following we provide, first, the solution in the case of 
a homogeneous and clean superconductor (where h = 0) and finally a solution for a ferromagnet/supcrconductor 
interface. 



1. Superconductor 

We now consider another limiting case: all potentials are set to zero apart from A^O. The Bogoliubov-De Gennes 
equation (|14|) again decouples into two sets of equations, one relative to spin f particles and spin [ holes 



A* -H J \ v{r) J \ v(r) 
and another, equivalent, relative to spin | particles and spin | holes. Considering plane wave solution of the form 



u if) \ _ ( V> 
v(f) 



e ik - 



(22) 



one finds the following dispersion relation: 



E = ±v/(fc 2 -^) 2 + |A| 2 , (23) 



which is plotted in Fig. [2] in the direction k z when |A| = fi/3 (solid line) and |A| = (gray line). We can notice that 
there is an energy gap in the energy spectrum for non-zero A as expected. It is useful to write down the solution for 
k z at a given energy E: 



k z = ±^±^j¥ = w , 



(24) 
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where fx = \i — k 2 — ky. This tells us that for an energy E above the gap there are four possible solutions denoted 
by q+ and q_ in Fig. [2] We shall denote excitations with wave- vectors k a as particle- like excitations and 

those with wave-vectors q a to be hole-like excitations. As we can see from Fig. [21 this choice means that particle-like 
excitations have their group velocity parallel to their momenta, whereas hole-like excitations have their group velocity 
anti-parallel to their momenta. This choice is made because the wave function associated with wave- vector k a (q a ) 
is predominantly particle- like (hole- like). This can be seen by substituting k a and q a in the amplitudes ip and </>, 
solutions of equation (f2Tj) : 

e iv/^ / Je 2 - |A| 2 




(26) 

where tp is the phase of the order parameter. Note also that since only the absolute values of the coherence factors ip 
and 4> are fixed {\ip\ 2 + \4>\ 2 = 1), there is an arbitrary choice in where to put the phase of the order parameter. 



2. Ferromagnet/ Superconductor (F/S) interface 

In this section we solve analytically the Bogoliubov-De Gennes equation for a F/S interface. We shall prove that the 
Andrccv reflection amplitude r a is suppressed when we are in the presence of an exchange field h. Before proceeding, 
note that, for a ballistic N/S interface, r a was first calculated by Blonder, Tinkham and Klapwijk (BTK) in Refi2£. 
r a was determined as a function of the energy of the quasi-particles and as a function of the strength of a barrier 
potential at the interface. The generalization of the BTK calculation to the case of a F/S interface was first derived 
in Ref^. 

For simplicity, let us consider a one-dimensional structure in which all potentials are set to zero apart from h = 
h 0(—x) z and A(x) = Ao 0(x), where h and Ao are constants. In this case, the Bogoliubov-De Gennes equation ^ 
can be decoupled into two equivalent equations, the first of which reads 

{-dl-n-h) A(x) \fuHx)\_^fuHx) 



A*(x) (d* + l x-h)J y V ±(x) J - E [vt(x) ) ■ (27) 

We now consider plane wave solutions of equation (|2"T)) and we solve the scattering problem. In the left-hand (x < 0) 
ferromagnetic region the wave function ipL^x) produced by a source of right-going spin j particles of unit flux at 
energy E can be written as 

u\x)\_ ( Jl+^o 



M*)=[ ) = [ Vl \ e ^ k ] » (28) 



r„ 



where k = y/E + JT+ h, q = \] —E + Jj, — h. Vk — 2fc and v q — 2q are the group velocities relative, respectively, to the 
particle wave-vector k and to the hole wave-vector q. r$ is the normal reflection amplitude, whereas r a is the Andreev 
reflection amplitude, which corresponds to the reflection of an incoming spin f particle into a spin J, hole. The Andreev 
reflection process^ consists of the coherent evolution of a particle-like excitation into a hole-like excitation. In the 
right-hand (x > 0) superconducting region the wave function ipR{ x ) is 



where ip^ and ± are the coherence factors, solutions of the Bogoliubov-de Gennes equation (|2Tj) for a superconductor 
(see sub-Section III C given by 

= ; 30 

E±y/E2-At 
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and |V' ± | 2 + l^l 2 — 1- to an d t a arc, respectively, normal and Andreev transmission amplitudes. In (]29[) we have 
defined: 



k + = \lv+\/ E ' 2 - A o > (31) 



= y / i- v /S 2 -A2, (32) 
dE _ 2k + (k\ -fi) 



dk z 



(33) 



and similarly for i> q+ . Since we are interested in the subgap solutions (i. e. when E < Ao), both fc + and q + arc 
complex and ipn(x) is an evanescent wave. The length £ over which the decay of \^>r(x) occurs can be defined by 
Re[ifc_|_£] = —1, so that £ = 1/Im[fc + ]. £ is known as the superconducting coherence length and for E = one has 
£ ~ It is therefore useful to define: k + = y/JT+Trj and q+ = y/ p — ir\ with 77 = yj Aq — E' 2 . The scattering problem 
is solved once ro, r a , to and t a are calculated. This can be done by matching the wave functions ipL and ipR and their 
derivatives at the interface (x = 0). We find 



2k(q + + k + ) 



(e + fc + )(0 - ?+) (J) - ( £ - + k+) (J) " 

In the case where H ~ and E 1 = 0, in the limit of small Ao//x, (|34[) reduces to 

/ , \ 1/2 

in such a way that the Andreev reflection probability R a = \r a \ 2 is a decreasing function of h. For small values of h, 
r a at E = can be approximated by 



(36) 



like in the N/S case. In real hard ferromagnets, like Co, the value of the exchange field, although high with respect 
to other ferromagnets, is one order of magnitude smaller than the Fermi energy, and therefore (|36p holds. In 3- 
dimensions, however, the picture changes since \i has to be replaced by \x — E ni where E n is the transverse kinetic 
energy relative to the n-th longitudinal mode. For quasi-particles approaching the interface at large angles, fi — E n 
can become comparable to or smaller than h and, as a result, the suppression of R a is enhanced with respect to the 
1-dimcnsional case. 

To conclude we remark, firstly, that r a contains the phase of the superconducting order parameter Lp. It is this 
which gives rise to interference phenomena for systems containing more than one superconducting interface with 
different phases. Secondly, Andreev reflection is the process by which electrical current can flow across a normal- 
metal/superconductor interface in the subgap regime. Since the change in momentum required for the Andreev 
process (Sp = h(k — q) ~ hkpE / 'p.) is much smaller than the change in momentum required for a normal reflection 
process (Sp = 2hk ~ 2Kkp), Andreev reflection is strongly favored at a clean normal-metal/superconductor interface 
for subgap energies. At the Fermi energy (E = 0) Andreev reflected hole retrace the trajectory of the incoming 
electrons (retroreflection) . During this process a charge of 2e is deposited into the superconducting condensate in 
the form of a Cooper pair, which would be carried away by a supercurrent. For energies greater that |A| normal 
transmission is dominant for clean normal-metal/superconductor interfaces. 



D. Scattering theory 

Electronic transport through a phase-coherent conductor can be studied in the same way as wave propagation into 
a scatterer. A conductor can be modeled as a scattering region connected to external leads which act as waveguides 
carrying a current of quasi-particles originating in external reservoirs. Scattering is assumed to be elastic and all 
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inelastic processes occur in such reservoirs where quasi-particles are distributed in energy according to the Fermi- 
Dirac function f(E) = [exp[(i? — /i) /kgT] + 1] _1 and characterized by a chemical potential /i = eV and a temperature 
T, e being the electronic charge. The conductance G of completely normal structures can be expressed in terms of 
scattering probabilitie o 27 ' 28 i 29 ' 30 yielding the Landauer-Biittiker formula: 

G = ^-Tr[tH] (37) 

for a 2-tcrminal system at zero temperature in the linear response regime. In Eq. (|37j) e 2 /h is the quantum of 
conductance, the factor 2 accounts for the twofold spin degeneracy and t is the transmission amplitude matrix for an 
electron to propagate through the one-dimensional conductor. 

For hybrid systems the scattering theory was first applied 3 -^ to a NIS interface (I being an insulating barrier of 
arbitrary strength) in Ref^. The total current Ins was calculated in the normal side of the junction in terms of the 
difference between incoming and outgoing quasi-particlc distribution functions yielding the so-called BTK formula: 

/oo 
[f(E - eV) - f(E)][l + A(E) - B(E)]dE , (38) 
-oo 

where A(E) and B(E) are the energy-dependent Andreev and normal reflection probabilities, respectively, and V is 
the bias voltage applied to the junction. In Eq. (|38|) iV(0) is the density of states at the Fermi energy, vf is the Fermi 
velocity, A is the effective cross-sectional area and / is the Fermi distribution function. A and B were determined by 
imposing matching conditions to the wavefunctions, solution of the Bogoliubov-de Gennes equation, at the N and S side 
assuming a (5-likc potential located at the interface. For low transmitting interfaces 7ns reproduces the conventional 
transfer Hamiltonian result for which the current is proportional to the density of states of the superconductor and 
to the transparency of the interface (Andreev processes arc completely neglected). Namely Ins is vanishingly small 
for sub-gap voltages and presents a pronounced peak at eV = A which rapidly decays asymptotically reaching the 
normal state current value. For highly-transmitting interfaces the I/V curve presents a large sub-gap current, due to 
Andreev reflection, which takes its maximum value when eV approaches A and thereafter decreasing to the normal 
state value. In addition it is also possible to determine that the length at which the quasi-particle evanescent wave 
penetrates into the superconductor is equal to Hvf/2A at energies close to the Fermi energy, i.e. of the same order of 
the BCS superconducting coherence length £ = hvp/irA. 

A more general treatment generalizing the Landauer-Biittiker formulae to superconducting systems was developed 
in Refsi 34 ' 35 ! 36 . It accounts for generic multi-terminal, multi-channel hybrid structures, even in the presence of 
ferromagnetism and spin-flip scattering. As mentioned before, in the absence of inelastic scattering (although a 
generalization to include this is available 3 ^) , dc transport is determined by the quantum mechanical scattering matrix 
S(E, H), which yields scattering properties at energy E, of a phase coherent structure described by a Hamiltonian H. 

Consider a scattering region connected to external leads which carry a current of quasi-particles originating in 
external reservoirs at chemical potential m = evi . In the presence of Andreev scattering, current-voltage relations for 
a phase-coherent scatterer connected to normal reservoirs were first derived in Refi^. If the structure is connected to 
external reservoirs by current-carrying leads with open scattering channels labeled by a set of quantum numbers n, 
then the S-matrix elements are defined through the relation: 

b n {E)=Y J S n ,n> (E,H) a n ,(E) (39) 

which connects the second quantization operators b n of quasi-particles leaving the scattering region through channel 
n to the operators a n < of quasi-particles entering the scattering region through channel v! . Both a n and b n satisfy 
anticommutation relations. The S-matrix satisfies the unitarity condition S^ 1 (E, TC) = S^E, Tt), due to quasi-particle 
conservation, and the time-reversibility condition S t (E,H) = S(E,TC*). In the presence of superconductivity it is 
useful to label the quasi-particle open channels in the leads using the following set of discrete quantum numbers: 
{i,a,aa}, where a = +1 characterizes particle-like excitations, a = —1 hole-like excitations, a is the open channel 



index in lead i and a = ±1 for spin up (down). In doing so, the quantity p" cr ' /3<T (E, Ti) = J2 a 



is the probability of reflection (if i = j) or transmission (if i =/= j) of a spin a' quasi-particle of type (3 in lead j to 

a spin a quasi-particle of type a in lead i. For a ^ (3 P^^ a is referred to as an Andreev scattering probability, 
while for a = /3, it is a normal scattering probability. In the presence of superconducting leads, we insist that all 
superconductors share a common condensate chemical potential fi. This is to avoid time-dependent order parameter 
phases varying at the Josephson frequency, which would render a time-independent scattering approach invalid. When 

the energy E is measured with respect to /is, the particle-hole symmetry Sy^fy &)(-^>^0 ~ a & ^{ia)'{j% {~E,H.) 
is satisfied too. 
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Because unitarity yields 



where z and j sum only over leads supporting open channels at energy E, we have 

i3,T',j 



and 



J2 P^' Pa (E,H)=N? a (E) 



(40) 



(41) 



(42) 



where N°" 7 (E) is the number of open channels for a-type a-spin quasi-particles of energy E in lead i. 
The total current due to all quasi-particles is the integral 



dEj2(a) 

CV(T 



°>P"'(E)f?(E) 



(43) 



where ffj(E) is the Fermi distribution of incoming quasi-particles of type j3 from reservoir j at temperature Tj and 
chemical potential fij : 



1 



(44) 



1 



and fi is the chemical potential of the condensate, when superconductors are present. At this point it is important to 
notice that all scattering coefficients can be computed in the presence of self-consistently determined order parameter 
A and normal potentials U a . In such a case one must take into account the effects on the system of changes in the 
applied potentials v/s. Such changes affect both the normal potential and the superconducting order parameter. For 
normal structures, self-consistency in the normal potential has been considered by Buttikcr et al. (Refi^) and it 
turned out that such effect can become important at large voltages. Calculations taking into account a self-consistent 
order parameter have been carried out for one dimensional structures in Ref . 39 i 40 ' 41 . These demonstrate that provided 



the currents are low enough with respect to the critical current, the matrix P i 
the application of finite voltages of order |A|. 
In the linear response limit, i.e. for (vj — v) 



0, we write down Ii in the form: 



h = E a ^ ( V J ~ v "> ■ 



can remain unchanged, even by 



(45) 



with 



7>(-i)?[^' (E){ 



P+°>+° { E)+P+°>-°(E) 



(46) 



At zero temperature, where — M; = 8(E), (|46|) reduces to 



*« ~ h 



(47) 



We conclude by noting that, at finite voltages, the chemical potential of the condensate /i should be determined 
self-consistently by imposing the conservation of the current J2i 1% = 0- From the expression of Ii in (|43|) it is clear 
that such a self-consistency condition involves integrals over all incident quasi-particle energies, which in turns require 
a knowledge of the S-matrix over a range of energies. In the case where superconductivity is present in one lead (say 
j) the solution is trivial, since \x coincides with one of the voltages vj. In the general case, however, the solution of 
the self-consistent equations is difficult, but for some simple structures may be determined by symmetry arguments. 
In a spatially-symmetric 2-probe system, for example, one can assume that = \(fX\ + 
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E. Two-probe differential conductance and conductance in more details 



We now apply the above formalism to derive two-terminal conductance formula?. We first consider the case of two 
non-superconducting leads (from now on, referred to as NN system) and then the case where one of the leads is a 
superconductor (NS system). Ii denotes the current entering the scatterer from the lead i 6 [1,2] and Vi is the bias 
voltage applied to it. Current conservation requires that 1% = —I\ = —I. 

In the NN system the conductance at zero-temperature and zero-voltage can be calculated by rewriting equation 
(1451) in matrix form as 



and considering its inverse 



I 

-I 



"it 

l 21 



'11 



'12 



l 21 "22 



"12 
1 22 



Vl — V 

v 2 - V 



(48) 



V2 



1 ( a 2 2 
d I -«21 



-ai2 
<<ii 



I 

-I 



(49) 



where a 
as G 



,1 



i 



a!; ; + ajj and d = 011022 — 012021 is the determinant. Since we are interested in the conductance defined 
, we want to get rid of the condensate chemical potential v. We can achieve this goal by subtracting the 



second equation in (|4"9"|) from the first one, obtaining 

G = - — 



d 



On + Oi2 + 2 1 + 022 

where the matrix elements , evaluated at the Fermi energy, are expressed as follows: 

,.2 



an = -r 

h 



7V T - i?J T - R l 1 



R 



1 1 



Rl l + R l n l 



(50) 



(51) 



a 2 i 



(Tr+T^-T^-T^) + (Tn+T^ 



T, 



1 1 



T, 



(52) 



and 



012 = -r 
h 



T: 



'i 1 



T, 



'II 



T, 



'1 1 



T, 



'I [ 



T, 



T: 



T, 



'1 1 



T, 



'1 1 



(53) 



022 



- R - R 



'IT 



R 



'TT 



R 



M l 



f N 'l 



- R - R 



R 



*\ 1 



R 



*\ 1 



(54) 



Here T CTCT 



= P^ CT = ^2i CT a ) IS t ne probability for normal (Andreev) transmission for an injected er'-spin 



quasi-particlc from the left-lead into a er-spin quasi-particle (hole) in the right-lead. For the primed quantities, the 
quasi-particles are injected from the right-lead into the left-lead. Note that T^ a and Rq 17 , for a ^ a' . and R a a a and 
T° a are non-zero only if spin-flip scattering is present. As a check of consistency it is easy to prove that in the absence 
of superconductors one recovers the well-known Landauer formula: 



G =~h 



T, 



1 1 



T, 



11 



T, 



i I 



T, 



: 1 



(55) 



When superconductors are present, it is interesting to note that their common condensate chemical potential \x can 
be determined imposing the conservation of the current which yields: 



(fl!2 + 022)^1 + («22 + «2l)^2 
811 + 012 + &21 + 022 



(56) 



In the NS system, we additionally insist that the superconducting lead and the superconductors in the scattering 
region share the same chemical potential: v 2 = v. Note that the relation I 2 = —I\ allows us to avoid the explicit 
calculation of the current in the superconducting lead (12), which will be a combination of quasi-particlc current and 
supcrcurrcnt. 
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The conductance at zero-temperature and zero- voltage is also easily derived. Equation ([48]) becomes trivial since, 
on the one hand, no transmission of quasi-particles to the superconductor is allowed at the Fermi energy (5i2 = 0) 
and, on the other, v-i — v — 0. Therefore (|48|) reduces to I = aw(v\ — v) which gives 

G = = j [(n} - i?J T - i?i T + RV + R l j) + (n± - Rl l - R l l + Rl l + R l a 1 )] (57) 

where the probability coefficients are calculated at the Fermi energy. 

It is very instructive to consider a structure composed of a superconducting scatterer of length Lg attached to 
two normal leads, one on the left and the other on the right. From Eq. P5|) one can derive the following 2-probc 
linear-regime zero-temperature conductance 4 ^: 



2e 2 

G = - 



T + T a 



2(R a R' & - T a Ti) 
i?a + RL +T a + TL 



(58) 



where T (T a ) is the normal (Andreev) transmission probability for quasi-particles injected from the left lead and 
arriving on the right lead. R (i? a ) is the normal (Andreev) reflection probability for quasi-particles injected from the 
left lead. Similarly T' & and R' a are Andreev scattering probabilities for quasi-particles injected from the right lead. For 
Lg 3> £ all transmission probabilities tends to zero, since the quasi-particles penetrate into the superconductor up to 
a depth of the order of £, so that the conductance reduces to the series of two interface resistances, namely l/2i? a 
and l/2i? a in units of h/(2e 2 ). A consequence of this is the fact that the overall conductance of the structure does 
not depend on Lg, but simply on the microscopic structure of the NS interfaces. In particular, for L$ — > oo although 
the resistance of the system attains an asymptotic finite value, the resistance per unit length (resistivity) vanishes, as 
it must be for a superconductor. 

According to Eq. the problem of determining the current- voltage characteristic is reduced to the calculation 

of scattering amplitudes. This can be done in different ways. In simple ballistic systems, for example, by solving the 
Bogoliubov-de Gennes equation piecewise in the different homogeneous regions composing the scatterer and imposing 
proper matching conditions to the wavefunctions at the boundaries between these regions. Complex structures such 
as disordered and diffusive conductors, heterostructures composed of different materials and multi-terminal systems 
can be dealt with too. Scattering amplitudes can be determined numerically, for example, by discretizing the system 
in real space within the tight-binding formalism. For diffusive wires and chaotic quantum dots a different approach 
can be successfully employed, namely the random-matrix theory (see, for example, the review papers of Refs . 43 ' 44 ). 
According to this theory, from the statistical properties of a class of matrices with random elements describing a certain 
physical system it is possible to extract the properties of the system. This can be applied to scattering matrices: the 
first problem consists in determining their random-matrix probability distribution (statistical ensemble). The second 
problem is to find the correlation functions of the transmission eigenvalues from which the transport properties can 
be derived. The conductance, for example, is simply given by 2e 2 /h /„ dTTp(T), where p(T) is the mean eigenvalue 
density. A complete solution to this problems has been found for diffusive wires and quantum dots. For hybrid systems 
consisting of a phase-coherent structure connected to superconducting leads, the random-matrix theory is based on a 
relationship which links the Andreev reflection to the transmission eigenvalues of the corresponding normal system 4 ^: 

4 e 2 T 2 

n— 1 v ' 

where T n is a transmission eigenvalue and N is the number of open channels 4 ^. From Eq. (|59p it is clear that hybrid 
systems can be studied making use of the results developed for normal conductors. 



F. Quasi-classical Green's function approach 



In this section we briefly underline another very important technique for dealing with hybrid systems. This approach 
applies when the characteristic length scale of the problem are large if compared with Fermi wavelength. It proves to 
be particularly useful for dirty structures where a sufficient amount of non-magnetic impurities is present so to make 
the motion of electrons isotropic. 

The quasi-classical theory of superconductivity is developed through quantum field theory methods and it is for- 
mulated in terms of Green's functions in the particle- hole space (known as Nambu space). Such Green's functions 
are denoted by a hat and, in addition to standard Green's functions G, they contain anomalous components F which 
describe the superconducting correlations: 

rtv r- r \-f G (^ r i' r 2) F(E,r u r 2 ) \ , v 

G^n.ra)-! fHe>F1)F2) G t ( £ jri , r2) ) • (60) 
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The starting point is the Gorkov equation^ for the Green's function of a bulk superconductor which was derived 
directly from the BCS Hamiltonian. The quasi-classical approximation consists in averaging the Green's function G 
over the rapid oscillations in the relative space coordinates p = rj — r2 and over impurities. As a results one obtains 
the quasi-classical Green's function g(E, r, n) which depends on the energy E, on the center-of-mass coordinate 
r = (ri + V2)/2 and the versor of the relative momentum n associated to the relative coordinate p. The Gorkov 
equation reduces to the following (Eilenberger) equation for g: 



v F d r g 

with 



-iEa 3 + A + ^-,g 
It 



, (61) 



( A* ) 



derived in Refs . 48 i 49 . Here vp is the Fermi velocity, A is the pair potential and r is the elastic scattering time. 
Furthermore the square brackets represent the commutator, while the angular brackets (• • • ) denote an averaging over 
the direction n. Note that Eq. (|61|) determines g up to a multiplicative constant and the following normalization 
condition must be applied: g 2 = 1. 

A further simplification can be realized in the presence of an isotropic impurity scattering potential which makes 
the motion of electrons diffusive (dirty limit). In this case g can be expanded in spherical harmonics and Eq. (|61|) 
reduces to the Usadel equation^ for the isotropic function g(E,r) = (g(E, r, n)): 

Dd r (gd r g) + [iEa 3 - A, g\ = (63) 

where D = v f t/3 is the diffusion constant and vp is the magnitude of the Fermi velocity. It is worthwhile to stress 
that both Eq. (|6"T|) and |63|) must be supplemented with a self-consistent condition for the pair potential A. The 
electrical current can then be calculated once the Green's function g is determined. 

It should be noted that the quasi-classical approximation does not allow to take into account nonuniformities 
that occur on the Fermi wavelength scale, such as boundaries, barriers and interfaces with other materials. It was 
shown, however, that this problem can be circumvented by applying proper boundary conditions to the quasi-classical 
Green's function. For the Eilenberger equation they were derived in Ref^i and for the Usadel equation in Refi^. It 
is interesting to note that such boundary conditions are obtained by making use of the connection between scattering 
amplitudes and Green's functions. For a more detailed treatment of the quasi-classical theory we refer the reader to 
a number of review papers on the subject (for the most recent see, for example, Refsi 10 i 53 ' 54 ' 55 ). 

An important progress in the quasi-classical theory was put forward by Nazarov who formulated it, within the 
dirty limit, in terms of a circuit theor y 56 i 57 i 58 . This constitutes a generalization of the classical Ohm's law to hybrid 
coherent nanostructurcs. These can be viewed as coherent networks consisting of nodes linked by connectors. The 
starting point consists in observing that when the second term in the Usadel equation (f6"3")) can be neglected, the same 
equation can be written in the form of a conservation law for a matrix current j (r)^: 

S r j(r) = 0; j(r) = aGd r G (64) 



where a is the normal state conductivity proportional to the diffusion constant D. The second equation in (|64[) 
resembles the local Ohm's law in a normal metal structure, while the first expresses the conservation of the current. 
The problem is finally defined by imposing additional boundary conditions, which sets the bias voltage established 
across the structure and describe interfaces in the conductor. By solving those equations the voltage distribution 
over the conductor and hence the local current density and total current can be determined. The ordinary circuit 
theory is formulated by discretizing Eqs. (|6"4")l plus boundary conditions through a finite clement approximation so 
that the structure can be separated into resistive elements pairwise connected. Current conservation in each node 
and current-voltage relations for each element set the voltage in each node allowing the determination of the total 
current. 

The extended circuit theory cannot be derived in such a simple way because of the matrix form of the equations 
and the fact that one deals with matrix currents and Green's functions. However the situation can be simplified due 
to the symmetry properties of the Green's function at zero energy so that Eqs. can be separated into two parts. 
The first one determines the equilibrium spectral local properties and the second one the non-equilibrium chemical 
potentials in each node which determines the electron propagation. These equations allows to formulate a circuit 
theory which can be conveniently condensed into a set of rules (see Ref^) for determining currents and potential 
differences of a coherent network of nodes and connectors. It is worth mentioning that: i) nodes are required to contain 
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enough disorder to prevent ballistic transmission between the two connectors attached to it. Furthermore the size of 
a node must be smaller than the coherence length, ii) connectors can be of three types, namely diffusive conductors, 
tunnel junctions or ballistic constrictions. The net effect of superconductivity is to induce a renormalization of the 
conductivity in the last two types of connectors. 

It is important to recall that the above circuit theory is valid when the second term in the Usadel equation can 
be neglected, i.e. for small temperatures and voltages 5 ^: kBT,eV <C A,hD/L 2 , L being the system size. The 
limitations described above are quite restrictive, namely the voltage and temperature dependence of the transport 
properties cannot be accessed. However, a circuit theory can still be formulated by discrctizing the Usadel equation 
(|63|l comprised of the second term. This can be done by discretizing the space into a set of nodes Xj such that 
the Green's function G(xi) = Gi in neighboring points are close to each other. A resistor is associated with each 
connection. The Eq. (|63[) can then be viewed as a conservation law in each node i (Y]u 1^ + Ii ca k a gc.i = 0) where now 
the total current is composed of Iik (which corresponds to the second expression in Eq. (|64p ) and an additional leakage 
contribution /leakage, i oc —ie 2 [Gi, H] associated to the second term of the Usadel equation ([63]) . The leakage current 
describes two processes, the first one, proportional to the energy e, is associated to the decoherence of electrons and 
holes due to wave- vector mismatch at finite energies ("leakage" of coherence). The second one, proportional to A, 
is responsible for the conversion between quasi-particles and Cooper pairs ( "leakage of quasi-particles" ) . In order to 
obtain an acc urate agreem ent with the full theory the distance between the nodes has to be smaller than the coherence 
length £ = \f DMax(A, e). It is important to notice that ballistic point contacts can be considered as a particular 
kind of connectors, but they need a specific treatment. To conclude, this new circuit theory can be formalized into a 
set of rules given in Refi 5 ^. 

G. Panorama on results in NS systems 

Superconducting hybrid nanostructures possess a rich physics related, in particular, to the interplay between co- 
herent transport of electrons and Andreev reflection. Physical systems investigated so far comprise NS junctions, 
mesoscopic SNS Josephson junctions, structures containing several superconducting islands and, more recently, ferro- 
magnet/superconductor (FS) junctions, where the interplay between spin-dependent transport and superconductivity 
can be investigated. The vast knowledge available in semiconductor technology has been also exploited by replacing in 
the above systems normal- metals with semiconductors (Sm). The possibility of using unconventional superconducting 
materials, such as d-wave superconductors and, more generally, high-T c superconductors, has been explored too. As 
a result, a variety of effects and phenomena have been reported so far, the most important of whose we shall review 
in the present section making use of the various theoretical approaches outlined in the previous sections. It is worth 
mentioning that these allow one to calculate, in addition to the electrical current, all transport properties of a given 
system such as thermal current, current noise and the whole full counting statistics of electronic transport. 

We shall start considering NS junctions, which present different properties depending on whether N is diffusive or 
ballistic and on the transparency of the barrier at the interface. Perhaps the first effect that received much attention is 
the zero bias anomaly (ZBA), which refers to a conductance peak observed at zero bias voltage in a low-transparency 
NS point contact 6 ^. This contrasts with the BTK result, where the conductance presents a minimum at zero bias, 
is particularly striking since the height of the peak is of the order of the normal state conductance. ZBA can be 
explained as an interplay between Andreev reflection at the interface with the superconductor and disorder scattering 
in the N side of the junction due to the presence of impurities. When the elastic mean free path in the N side of the 
junction is smaller than the junction size, electrons have the chance to be scattered back to the S interface many times, 
finally undergoing an Andreev reflection. The net result is that Andreev reflection processes occur with a much larger 
probability with respect to the case of absence of impurities. In other words the low bias conductance is determined 
by complex interference effects which produce an enhancement by several orders of magnitude. Such an effect dies 
away for energies larger than the Thouless energy which sets the scale for particle-hole dephasing. Although ZBA 
was first understood within the quasi-classical approach, it was then confirmed using scattering methods and also in 
the tunneling Hamiltonian approach. This effect is also known as reflectionless tunneling, since it can be explained in 
terms of the disorder-induced opening of tunneling channels. 

Non-monotonic behaviors of the conductance as a function of voltage and temperature in NS junctions with clean 
interfaces were also reported. When a diffusive wire is connected to a superconducting reservoir through a highly 
transparent contact the zero-temperature, zero- voltage average conductance equals the normal-state conductanc e 45 ' 61 . 
Even though this can be proved rigorously, intuitively it is the consequence of two facts, namely that the Andreev 
reflection effectively doubles the length of the normal metal conductor and that the conductance in a NS point contact 
is doubled with respect to the normal-state one. However it was 

foundSLSSiSaiS^S^ that a conductance peak appears 
either for voltages or temperatures of the order of the Thouless energy HD/L 2 . This phenomenon, known as reentrance 
effect, was experimentally proved in Refsi 67 ' 68 ' 69 . 



14 



Another interesting phenomenon is the Andreev intcrfcromctry which is realized in a hybrid structure containing 
at least two superconducting islands. This is based on the fact that when a quasi-particle Andreev reflects from a 
NS interface, the phase of the outgoing excitation is shifted by the phase of the superconducting order parameter. 
In the presence of two superconducting islands with order parameter phases <f>\ and </>2, the transport properties are 
oscillatory functions of the difference <p\ — <p2 ■ The conductance of individual samples was found to be 27r-periodic in 
the phase difference. 

During the last few years the interplay between Andreev scattering and the ferromagnetic order has been addressed 
in numerous studies of electronic transport properties in nanostructurcs containing both fcrromagnets and supercon- 
ductors. The interest comes from the fact that the electron and the hole involved in the Andreev process must belong 
to opposite spin band. This produces a suppression of the current flowing through a FS interface, which depends on 
the value of the exchange field characteristic of the ferromagnet. This phenomenon has been used for estimating the 
spin-polarization of different kind of ferromagnet a 70 ' 71 . 



III. ELECTRONIC COOLING IN SUPERCONDUCTING NANOSTRUCTURES 



A. SIN refrigerators 



An important research topic is represented by heat transport in superconducting nanostructurcs. In particular, heat 
transport through SN interfaces can be successfully applied to microcoolin g 72 ' 73 . In order to gain some insight into 
this problem it is instructive to consider the simplified energy band diagram of a NIS tunnel junction (I is an insulating 
barrier) biased at voltage V, as schematically depicted in Fig. [3fa). The physical mechanism underlying electronic 
cooling is rather simple. As discussed in section lll C| when a normal metal is brought in contact with a superconductor, 
quasiparticle transport is effective only at energies larger than the superconducting gap A. As a consequence of this 
"selective" tunneling of hot particles, the electron distribution function in the N region becomes sharper, i.e., the 
electron effective temperature in N is lowered (even in the regime when electrons are thermally decoupled from the 
lattice): the SIN junction thus behaves as an electron cooler. This situation can be experimentally accessed when 
transport is dominated by quasiparticle dynamics, i.e., in NIS tunnel junctions. 

To highlight the impact of junction transmissivity in governing heat flux (cooling power) jq^ through a voltage 
biased SN contact (see Fig. [3Jb)) let us consider Fig. EJc), where the maximum cooling-power surface density jqa 
versus interface transparency T> is calculated at different temperatures, jqa is a non monotonic function of interface 
transmissivity, going to zero both at low and high values of T>. In the low transparency regime it turns out to be 
linear in 23Zi, i.e., the electron transport is dominated by single particle tunneling. However, by increasing interface 
transmissivity, two-electron Andreev tunneling^ begins to dominate quasiparticle transfer across the junction leading 
to enhanced conductivity but suppressing the heat flow extraction from the N metal, jqa is maximized at an optimal 
T> value, which in addition depends on the temperature. In real SIN structures for cooling applications (i.e., those 
exploiting tunnel junctions), V is typically in the range 10~ 6 -j- 10~ 472 i 73 corresponding to junction specific resistances 
(i.e., the product of the junction normal state resistance and the contact area) from tens to several thousands fl^m 2 . 
This limits the achievable jqa to values much lower than those expected from theory, in particular to some pW//jm 2 . 
Barrier optimization seems to be still nowadays a challenging task. 

In the tunnel limit (that is particularly relevant in the experiments), the cooling power is given by 

3Q(Y) = I dEEN s (E) [f N (E)-f s (E)} , (65) 

where E = E — eV , fs(N)(E>) = fo(E, Tg^) is the Fermi-Dirac distribution of quasiparticles at temperature T$ in the 
superconducting reservoir (and of electrons at temperature T e in the N electrode) , Rt is the normal-state resistance of 
the junction and Ns(E) = \Re[E/^E 2 — A 2 ]| is the BCS normalized density of states of the superconductor. Figure 
E](d) shows jq calculated from Eq. (|6"5")) versus bias voltage at T = T e = T$ = 0.3 A/Zcg. The cooling power is a 
symmetric function of the bias voltage V, i.e., Jq{V, T) = ]q{— V, T), and, at this temperature, it is maximized around 
V w ±0.8 A/e, as indicated by red-dashed arrows in the figure. Note that when jq is positive, it implies heat removal 
from the N electrode. Furthermore, for each temperature, there is an optimized voltage that maximizes jq. The inset 
of Fig. 03d) displays the heat current versus temperature calculated at each optimal bias voltage. The cooling power 
is maximized at T w 0.3 A/ks = T opt (indicated by the blue-dashed arrow) where it reaches jq ~ 6 x 10~ 2 A 2 /e 2 i?T, 
decreasing both at lower and higher temperatures. 

The final T e is determined from the balance among several factors that tend to drive power into the electron 
system. The majority of SIN microrefrigerators realized so far operate in a regime where strong inelastic electron- 
electron interaction drives the system to (gMasi)cquilibrium, where the electron system can be described by a Fermi 
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function at a temperature T e that differs, in general, from that of the lattice (T p /j). At the temperatures of interest 
(i.e., typically below IK), the predominant contribution comes from electron-phonon scattering that exchanges 
energy between electrons and the lattice phonons. This rate of exchanged energy is given b y 75 1 76 

f Q - ph = EV(T^ h -T!), (66) 

where £ is a material-dependent parameter of the order 10 -9 WK -5 /xm -3 and V is the volume on the normal electrode. 
The refrigerator cooling power {jq e ^ r ) can be defined as the maximum power load that the SIN device can sustain 
while keeping the N region at temperature T e , and can be generally expressed as 

3q fr = 3Q ( V ' T e, T ph ) - f Q ~ ph (T e ,T ph ). (67) 

The minimum achieved temperature T e thus fulfills the condition 2q^ T = 0- 

The first observation of heat extraction from a normal metal dates back to 1994Z 3 -, where cooling of electrons in Cu 
below the lattice temperature was demonstrated using an AI/AI2O3/CU tunnel junction. An important improvement 
was made two years later, still in the AI/AI2O3/CU system, by Leivo et al&, which recognized that using two SIN 
junctions in series arranged in a symmetric configuration (i.e., in a SINIS fashion) would have led to a much stronger 
cooling effect. In this configuration a reduction of the electron temperature from 300 mK to about 100 mK (i.e., close 
to the achievable optimum) was obtained. Later on, several other experimental evidences of electron cooling in SINIS 
metallic structures were reporte d 77 ' 78 ' 79 ' 80 ' 81 ' 82 ' 83 ' 84 ' 85 ' 86 ' 87 ' 88 ' 89 . In these experiments SIN junctions are used to alter 
the electron temperature in the N region as well as to measure it. In order to measure the temperature, the N region 
is normally connected to other SIN contacts (i.e., "probing" junctions) that act as thermometers, suitably calibrated 
by varying the bath temperature of the cryostat. 



B. SF refrigerators 

As discussed in Sec. IIII A[ decreasing the contact resistance Rt is not a viable route to enhance heat current in 
SIN junctions. As a matter of fact, the increase of the Andreev reflection contribution to the electric transport^ 
(see Fig. Ufa)) strongly affects and limits thermal current through the system. A promising scheme to surmount this 
problem is to exploit, instead of an insulating barrier, a thin ferromagnetic (F) layer in good electric contact with 
£3. As already discussed in section III C 2| depending on the degree of spin polarization V of the F layer, strong 
suppression of the Andreev current may occur at the SF interface^ 6 .. In the limit of large V, the subgap current is 
thus strongly suppressed while efficient hot-carrier transfer leads to a considerable heat current across the system. 
The authors of Refi^ - focused on a ballistic three-dimensional SF junction (as shown in Fig. (He)) where the F region 
is described by the usual Stoner model^ 6 .; furthermore, the analysis was carried out along the lines of a voltage biased 
SN junctio n 25 ' 74 but generalizing heat transport to a spin-dependent superconducting structure. In the following we 
shall give the main results of such analysis. 

FigureUfd) shows T e as a function of bias voltage for two starting lattice temperatures (i.e., the system temperatures 
at V = 0) Tph = 200 mK and T p h = 300 mK. The calculation was performed for half-metallic (V = 1) SF (filled-green 
dots) and for SIN (filled-blues triangles) microcoolers, assuming Al as superconductor (Aai = 180 /j,eV). In particular, 
here it was assumed a junction area A = 0.1 /-tm 2 , an active volume of the FN side V = 0.5 /1m 3 , S = 2 x 10~ 9 WK 
" 5 /jm" 3 (typical of Cu), and a SF normal-state resistance R^ s = 28 mfl (assuming a Fermi energy £p = 5 eV), 
and for the SIN-junction specific contact resistance a value 0.4 Ml /mi 2 , corresponding to high-quality CU/AI2O3/AI 
junction o 72 ' 73 . Figure [3|d) shows the remarkable T e reduction provided by the FS cooler with respect to the NIS 
cooler. Even starting from 200 mK, for which the tunnel junction provides its largest cooling effect (a temperature 
reduction of about 10% at eV ~ A), the SF cooler yields T e of the order of 10 mK (a temperature reduction of about 
95%). This marked difference stems from the high contact resistance of the SIN junction that strongly affects its 
performance and must be compared to specific contact resistances as low as 10 -3 f2/itm 2 that are currently achieved 
in highly transmissive SF junctions 7 ^. 

The junction specific cooling power is shown in Fig. EJa), where jqa (evaluated at each optimal bias voltage) is 
plotted versus V at different bath temperatures. Notably, for a half- metallic ferromagnet at T = 0.4 A/fc^ cooling 
power densities as high as 600 nW//im 2 can be achieved, i.e., about a factor of 30 larger than those achievable in 
SIN junctions at the optimized interface transmissivity (i.e., T> ~ 3 x 10~ 2 at T ~ 0.3A/fcs, see Fig. Etc)). The 
impact of partial spin polarization (V < 1) in the F region is shown in Fig. [5jb) where the normalized heat current 
jq versus bias voltage across the junction is plotted for some V values at T = 0.4 A/fcs. For each V value there 
exist an optimum bias voltage which maximizes jq. In the limit of a half-metal ferromagnet (i.e., V = l) 92 i this 
value is around V op t — A/e. In the inset of Fig. El (b) is displayed jq calculated at each optimized bias voltage 
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against temperature for various V values. The heat current is maximized around T = T opt ~ 0.4 A/ks, rapidly 
decreasing both at higher and lower temperatures The inset of Fig. [5ja) shows the junction coefficient of performance 
(77) calculated at the optimal bias voltage versus temperature for various V values. Notably, for V = 1, 77 reaches 
~ 23% around T = 0.3 A/fcs and exceeds 10% for V = 0.96. In addition, in light of a possible exploitation of this 
structure with a N region (i.e., a SFN cooler) it turned out that cooling effects comparable to the SF case can be 
achieved with F- layer thickness of a few nm (i.e., of the order of the magnetic lenght). 

The above given results point out the necessity of strongly spin-polarized ferromagnets for a proper operation of the 
SF cooler. Among the predicted half-metallic candidates it is possible to indicate Cr0 2 94 i 95 ' 96 , for which values of V 
in the range 85 ~ 100% have been reported 9 ^^^, (Co^Fe^S^, NiMnSb^, Sr 2 FeMoOe±£3 an d NiMnV^. 
So far no experimental realizations of SF structures for cooling applications were reported. The main technological 
problems seem, still nowadays, to avoid degradation of ferromagnetic properties at the SF interface and a proper 
integration with suitable superconductors. 



IV. COULOMB BLOCKADE IN HYBRID STRUCTURES 



A. Coulomb blockade and superconductivity 



The astonishing progress in modern technology has made it possible to fabricate in a controlled way tunnel junctions 
with capacitances of the order of C = 10~ 15 F and below. In this case even the charging energy associated with a 
single-electron accumulation at the junction, Eq = e 2 /2C, can be of the order of several Kelvins. This implies that 
the properties in the sub-Kelvin regime are strongly affected by the presence of charging effects. Indeed charging 
energy allows the control of the electron number of small islands with precision e or 2e in the case of superconducting 
devices. Adding one electron to a small superconducting island puts it into an excited state with an energy exceeding 
the gap. Only when a second electron is added, can both recombine to form a Cooper pair. This leads to parity 
effects, i.e. to different properties of the superconducting systems depending on the number of electrons being even 
or odd. 

The relevant mechanisms for charge transport in superconducting single-electron devices are quasi-particle tunnel- 
ing, two-particles tunneling and coherent tunneling of Cooper pairs. All of them are strongly affected by charging. 
Tunneling of quasi-particlcs is very similar to the case of normal metals but taking into account the modified BCS 
density of states. Two-particle tunneling is a higher order process related to the Andreev reflection discussed previ- 
ously. Both two-particle tunneling and coherent oscillations of Cooper pairs dominate transport at very low voltages. 
Moreover, the interplay of charging and coherent tunneling has a deep significance because it leads to the possibility 
of observing macroscopic quantum dynamics in these systems. The charge and the phase difference in a Josephson 
junction, although being macroscopic degrees of freedom, are quantum mechanical conjugate variables. Therefore the 
eigenstates are in general superpositions of different charge states. In addition the combination of coherent Cooper 
pair tunneling and quasi-particle tunneling, leads to a variety of structures in the I-V characteristic. 

There are numerous examples of systems where single-electron properties and superconductivity have been studied 
theoretically and realized experimentally. A prototype example which we will analyze is the NSN-SET transistor 
where a superconducting island is connected to two normal conducting leads by means of tunnel junctions. This 
example allows to discuss both quasi-particle and two-particles tunneling. We then discuss some basic properties of 
the coherent oscillations of Cooper pairs which will be important for the discussion on the implementation of quantum 
computers by means of superconducting nanocircuits. 

As prototype examples to discuss charging effects in superconducting nanostructures we consider the single-electron 
(SET) transistor shown in Fig. [5^ and the single-electron box of Fig.O). Here we concentrate on the transistor, while 
the box will be analyzed in the section devoted to quantum computation. The charging energy of the SET transistor 
depends on the electron number in the central island and on the applied voltages. The central island is coupled via two 
tunnel junctions to a transport voltage source, V = Vl — Vr, so that a current can flow. The island is, furthermore, 
coupled capacitively to a gate voltage Vq. The charging energy of the system depends on the integer number of excess 
electrons n = ±1,±2, ... on the island and on the continuously varied voltages. Elementary electrostatics yields the 
"charging energy" 

E ch (n,Q G ) = ^§^. (68) 

Here C = + Cr + Cq is the total capacitance of the island. The effect of the voltage sources is contained in the 
gate-charge Qg = CqVq + ClVl + CrVr. Similar expressions hold for the single-electron box. 

In a tunneling process which changes the number of excess electrons in the island from n to n + 1, the charging 
energy changes. Tunneling from the left junction to the island is possible at low temperatures only if the energy in the 
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left lead, eVi,, is high enough to compensate for the increase in charging energy eVL > E c h(n + 1, Qq) — E c h(n, Qg)- 
Similarly, tunneling from the island (transition from n+ 1 to n) to the right lead is possible at low temperatures only if 
E c h(n+1, Qg) — E c h(n, Qg) > cVr. Both conditions have to be satisfied simultaneously for a current to flow through 
the transistor. If this is not the case, the current is exponentially suppressed realizing the so-called Coulomb blockade. 
Varying the gate voltage produces Coulomb oscillations, i.e. an e-periodic dependence of the conductance on Qg- 
To further understand the characteristic of a SET transistor we need to determine the tunneling rates associated to 
electron tunneling. An electron tunnels from one of the states k in the left lead into one of the available states q in 
the island, thereby changing the electron number from n to n + 1, with rate tli- Such a rate, which is calculated by 
means of the Fermi golden rule, can be expressed in a transparent way as 



1 1 

-I qp (0E ch e) — 

e cxp(dE; ch /fcBT) - 1 



7li (n, Qg) = -I qP (SE ch /e) /X1P /u _^ — - . (69) 



The function / qp (V) is the quasi-particle tunneling characteristic (see e.g. Refji), which is suppressed at voltages below 
the superconducting gap(s) and 5E c h is the charging energy difference between the states with n and n + 1 electrons 
in the island. Charging effects reduce the quasi-particle tunneling further. At low temperatures, k-gT <C |<5-E c h|, 
a tunneling process which would increase the charging energy is suppressed. This phenomenon is called Coulomb 
Blockade of electron tunneling. At zero temperature the rate is nonzero only if the gain in charging energy compensates 
the energy needed to create the excitations <J£7 c h + 2A < 0. 

The rates describe the stochastic time evolution of the charge of the junction system. For the theoretical analysis 
a master equation approach can be used. Several examples of the current- voltage characteristics of normal metal can 
be found in Refs.— an important characteristic is the e-periodic dependence of the current and conductance on the 
applied gate charge. 

As dictated by Eq. (|69[) . the rate for quasi-particle tunneling are exponentially suppressed below the supercon- 
ducting gap. In this case higher-order processes involving multi-electron tunneling can play a role. Indeed, for NS 
interfaces, there still exists a 2-electron tunneling process, in many respects similar to the Andreev reflection discussed 
in the previous section and for this reason denoted as Andreev tunneling. The rate for Andreev tunneling taking into 
account charging effects was discussed in RefJ^. Andreev tunneling is a second-order coherent process. In the first 
part of the transition one electron is transferred from an initial state, e.g. k | of the normal lead, into an intermediate 
excited state q | of the superconducting island. In the second part of the coherent transition, an electron tunnels 
from k' I into the partner state —q J, of the first electron, such that both form a Cooper pair. The final state contains 
two excitations in the normal lead and an extra Cooper pair in the superconducting island, and the rates reads: 

7liKQg) = -T^ JTrTf\ 7- (70) 

Ae z cxp(dE ch ,2/kBT) - 1 

Note that the functional dependence of this rate coincides with that for single-electron tunneling in a normal junction, 
Eq. (p39"|) . with a linear I — V characteristics. Hence Andreev reflection is subject to Coulomb blockade like normal- 
state single-electron tunneling^ with the exception that the charge transferred in an Andreev reflection is 2e, and 
the charging energy changes accordingly. The effective Andreev conductance is of second-order in the tunneling 
conductance and, as shown in Refp^Si, is sensitive to spatial correlations in the normal metal, which can be expressed 
by the Coopcron propagator. 

In a normal- metal electron box (see Fig. \^jp), by sweeping the applied gate voltage the the number of electrons on 
the island in unit of e and, as a consequence, the voltage of the island shows a periodic saw-tooth behavior-^. The 
periodicity in the gate charge Qg is e. If the island is superconducting and the gap A smaller than the charging energy 
Eq, the charge and the voltage show at low temperatures a characteristic long-short cyclic, 2e-periodic dependence 
on Qq- This effect arises because single-electron tunneling from the ground state, where all electrons near the Fermi 
surface of the superconducting island are paired, leads to a state where one extra electron - the "odd" one - is in an 
excited stated. In a small island, as long as charging effects prevent further tunneling, the odd electron does not 
find another excitation for recombination. Hence the energy of this state stays above that of the equivalent normal 
system by the gap energy. Only at larger gate voltages can another electron enter the island, and the system can 
relax to the ground state. This scenario repeats with periodicity 2e in Qg, as displayed in Fig. [71 

Parity effects appear at temperatures lower than a crossover temperature T c g which, in the experiments, is typically 
of the order of 10 — 30% of T c . At low temperatures the even-odd asymmetry has been observed in the single- 
electron bosii^ as well as in the I-V characteristics of superconducting SET transisto r m i 112 i 113 . However, at higher 
temperatures, above a cross-over value T cr <C A, the e-periodic behavior typical for normal-metal systems is recovered. 



18 



B. NSN Transistor 

The analysis sketched above can be used to analyze the I-V characteristics of SET transistors with a superconducting 
island, the so-called NSN transistor. The interplay of parity effects and Andreev tunneling makes the current-voltage 
characteristics quite rich. If the energy gap is smaller than the charging energy the important processes are single- 
electron tunneling processes in the left and right junction, causing transitions between even and odd states. In the 
opposite limit of a superconducting gap bigger than the charging energy, the odd states have a large energy and 
Andreev tunneling, with rate given by Eq. (|70j) . becomes importanbaJai. 

At low temperatures, and superconducting gap larger than the charging energy, a set of parabolic current peaks is 
found centered around the degeneracy points Qq = ±e, ±3e, . . ii2£ 

l\8Q Gl V)=G^(v-^Y( V -^) ■ (71) 

Here SQq = Qq — e for Qq close to e, and similarly near the other degeneracy points. At larger transport voltages, 
single-electron tunneling sets in, even if A > Eq , and Andreev reflection gets "poisoned"^££. This occurs for 

V > V poison = 1(ec-^+A^ . (72) 

Fig. [5] shows the current- voltage characteristic of a NSN transistor with A > EqP^. At small transport voltages 
the 2e-pcriodic peaks due to Andreev reflection dominate; they get poisoned above a threshold voltage. The peaks at 
larger transport voltages arise from a combination of single-electron tunneling and Andreev reflection. The shape and 
size of the even-even Andreev peaks and some of the single-electron tunneling features at higher transport voltages 
agree well with those observed in the experiments of Hergenrother et alpii^. 

V. CONCLUSIONS 

In this review we tried to present an overall description of mcsoscopic hybrid structures. The field is so vast that 
the choice of the topics was surely biased by our personal taste and by our field of investigations. To summarize, 
we have started this review by illustrating the transport properties of hybrid systems formed by putting into contact 
superconductors with normal metals and fcrromagncts. In Section II we have introduced the Bogoliubov-dc Gcnncs 
equation in order to discuss transport properties in terms of Andreev processes. We have then provided a derivation 
of the Bogoliubov-de Gennes equation starting from the BCS theory and given the solutions of such an equation for 
two different systems, namely a homogeneous superconductor and a ferromagnet/ superconductor junction. After 
that we have introduced the two most important theoretical approaches to the study of transport properties in hybrid 
systems. We chose to give a more detailed treatment of the scattering theory, deriving general conductance formulas in 
terms of scattering matrices and giving explicit expressions for the case of two-probe systems. The quasi-classical 
Green's function approach has also been discussed, although in less details, and a review of the most important results 
has been provided. In Section II we have considered the heat transport in SN and SF structures to discuss the very 
efficient electronic cooling properties of hybrid systems. Section III has been devoted to the Coulomb blockade, which 
is present in small metallic grains where charging effects are important, in hybrid systems. We have discussed parity 
effects in single-electron transistors containing a superconducting island and their consequences on the quasi-particlc 
tunneling rates. I-V characteristics has been furthermore illustrated. 

There are several important topics which were left out, the interested reader may find additional references in the 
reviews quoted in the Introduction. 

It is quite hard to draw conclusions and even harder to advance ideas about possible future directions of the 
field. Nevertheless we take this risk and briefly outline what we believe (remaining conservative) will be the future 
development of the field. A direction which is not fully studied up to now is the role of superconductivity in molecular 
electronics. Transport through molecules using superconducting electrodes may give additional insight due to the fact 
that one can perform both single-particle and Cooper pair spectroscopy. 
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FIG. 1: a) Energy spectrum of quasi-particle excitations in a superconductor: quasi-particle energies E are plotted as a 
function of their wave-vectors k. Superconducting correlations give rise to a gap of amplitude |A| above the Fermi level set 
by the horizontal axis. Typical values of A go from hundreds of fieV to a few meV. b) In a normal reflection process, which 
occurs for example at an interface with an insulator (I), an electron is scattered back into the normal metal conserving the 
component of the momentum parallel to the interface, c) In an Andreev reflection, occurring at a NS interface, the incoming 
electron is reflected back as a hole. In this case the component of the momentum parallel to the interface is reversed, so that, 
at the Fermi energy, the reflected hole retrace back the trajectory of the impinging electron. 
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FIG. 2: Dispersion curve relative to a uniform superconductor for A = 1/3/^ (solid line) and A = (grey line), taking A = 1, 
jj, — 3 and k x — k y — 2n /20. For A / an energy gap of amplitude |A| appears in the spectrum. 
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FIG. 3: (a) Energy band diagram of a biased superconducting tunnel junction. More energetic electrons (e) can more easily 
tunnel into the superconductor. The electron gas in the N electrode is thus cooled, (b) Schematic representation of a voltage 
biased SN junction. Upon biasing the structure, a heat current jq flows across the system, (c) Maximum specific cooling power 
Jqa versus interface coefficient of transmission T> at three different temperatures, (d) Dimensionless cooling power jq of a SIN 
contact versus bias voltage V at T — 0.3 A/fes- The inset shows jq calculated at the optimal bias voltage versus temperature. 
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FIG. 4: (a) Schematic description of Andreev reflection at a SN contact, (b) Schematic representation of the principle of 
operation of a SFN electron cooler. For V = 1, Andreev reflection is hindered by the absence of available states for reflected 
holes (h). This subgap electron-transport suppression mechanism allows the operation of the microrefrigerator in the presence 
of efficient carrier transfer to S. (c) Sketch of the SFN electron cooler. The FN contact is supposed to be a highly transmissive 
electric contact, (d) Electron temperature T e vs bias voltage V for a half-metallic (V = 1) SF and SIN microcoolers for two 
starting bath temperatures T p h at V = 0. Adapted from Ref.— and Ref.— . 
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FIG. 5: (a) Maximum cooling-power surface density jqa versus spin polarization V for various temperatures. The inset shows 
the coefficient of performance r\ calculated at the optimal bias voltage versus temperature for some V values, (b) Dimensionless 
cooling power jq vs bias voltage at T = 0.4 A/fcs for some V values. The inset shows the heat current calculated at the optimal 
bias voltage against temperature for several V values. Adapted from Ref.— and Ref.— . 
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FIG. 6: a) The SET transistor, b) The Cooper-pair box. 
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FIG. 7: The charging energy of a superconducting single-electron box as a function of the gate voltage shows a difference 
between even and odd numbers n of electron charges on the island. Accordingly the average island charge (n) is found in a 
broader range of gate voltages in the even state than in the odd state. 




